^ = at a Sharp Semiconductor/Insulator Interface. 

Is This Correct? 



Leonid S. Braginsky* 
Institute of Semiconductor Physics, 630090, Novosibirsk, Russia 

(February 1, 2008) 



Abstract 

The generalized boundary conditions for the envelope wave function that take 
into account the real structure of an interface were used to investigate the 
hole spectrum of the semiconductor quantum dot embedded in an insulator 
matrix. An essential influence of the interface levels, which could exist at the 
top of the valence band, on the hole spectrum has been demonstrated. It is 
found that usually applied boundary conditions, wherein all components of 
the envelope wave function vanish at the interface, can be used only in the 
absence of the interface levels close to the band edge. 
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To determine the electron states in semiconductor quantum dots, the envelope function 
approximation is applied.0 The envelope wave functions are usually supposed to be vanished 
at the interface in order to describe the quantum confinement. If the intervalley or interband 
degeneracy of the electron spectrum occurs in the semiconductor, then each component of 
the envelope wave function is assumed to be vanished at the interface. The last statement, 
however, needs to be justified. 

Indeed, the boundary condition = for the proper (not envelope) wave function 
arises from the solution of the Schrodinger equation for the infinitely high step-like potential 
barrier. Such a barrier can't be considered in the framework of the envelope function ap- 
proximation wherein the potential must be smooth on the scale of the lattice constant. This 
means that the boundary condition \I' = is justified if the potential that restricts the elec- 
tron movement is smooth on the scale of the lattice constant, but sharp on the scale of the 
electron wavelength. Perhaps, this happens at the contacts of chemically similar materials 
(e.g., GaAs/AlAs),0 but not at the contacts of too different materials (e.g., semiconductor 
microcrystals embedded in the glass matrix). 

The interface influence on the electrons in the quantum dot has to be more complicated 
if the simple boundary condition = doesn't hold. It is well known that band mixing 
exists at the interface in the heteroj unctions of AniBy semiconductors.ili However, this is 
impossible at the plane interface where = 0. The bound electron states are also possible 
at the interface owing to the interband mixing;0 in these states the electron wave function 
behaves roughly as ip oc [exp (— 7ir) — exp (—72?")], i.e., the wave function vanishes at the 
interface and far from it, but has an extremum at a certain distance from the interface. 

The band mixing arises when the crystal symmetry that leads to the interband degener- 
acy in the bulk disappears at the interface. The lattice constant is the size that is charac- 
teristic for such symmetry. Therefore the proper boundary conditions for the envelope wave 
functions should take into account the real structure of the interface. 

In this paper we propose simplest boundary conditions that take into account these 
subtle details of the interface influence. We find the conditions under that the boundary 
condition = is applicable at a sharp interface. 

Let z = he the plane interface between a semiconductor {z > 0) and an insulator 
{z < 0). Assume the two-fold degeneracy for the electron band of the semiconductor and 
nondegenerate electron band in the insulator. Then the boundary conditions for the envelope 
wave functions of the electrons in these bands can be written as follows:i 

^i(r°) = bn^iini), 

*2(r2°) = fe22^K^22), (1) 
&31^i(t-3i) + &32^2(l-32) = ^K^s)- 

Where \E'i and \I'2 are the envelopes that relevant to the degenerate band of the semiconduc- 
tor, and \E'; is the electron envelope in the insulator. To obtain the parameters bij and Tij of 
the boundary conditions (|1|), the Shrodinger equation has to be solved in the narrow (about 
a few lattice constants) region at the interface. It is impossible at an arbitrary and rather 
imperfect interface. Nevertheless, these parameters are independent of the electron energy; 
they characterize the interface, and estimations of their values (the small width of the in- 
terface region is the fact that is important for these estimations) are bij ~ 1 and \Tij\ ~ a, 



2 



where a is the lattice constant. Thus, the boundary conditions (|I]) take into account the 
real structure of the interface. 

We assume the effective-mass approximation holds in the bulk of each material, so that 
\l/(r) = \l/(0)+r\I''(0). The large bands offset at the interface restricts the electron movement. 
If so, then oc exp (7/2;) and = •yi'^i, where the 7/ value can be considered as independent 
of the electron energy. Eliminating from the Eqs. (|lD yields 

^i(rii)+6i2^2(ri2) =0, (2) 
^21*1(^21) + ^2(r22) =0, 

where bij ~ bij and fjj ~ Xjj are known functions of bij, Tij, and 7^. To ensure the probability 
flux conservation at the interface, we have to assume 

^12(7^22 - ^2) _ hiihi - T21) , s 

mi m2 ' 

where mi and m2 are effective masses of the appropriate bands. 

The Eqs. (H) are the general form of the boundary conditions that should be written 
instead of \l/ = at a sharp semiconductor/insulator interface. The most general boundary 
conditions that are applicable at such interface have been considered in Ref.l3. Ours, Eqs. (^, 
hold in the effective-mass approximation. This approximation has been used ini to obtain 
the Eqs.(0) and to estimate the parameters bij and Ttj. 

It is important that the boundary conditions (|]) are nonlocal; they relate the envelopes 
at the different points fij near the interface. However, the mean width of the "nonlocality 
region" is small in comparison with the electron wavelength A (|fjj| ~ a <^ A) . To under- 
stand consequences of this nonlocality, let, at first, assume fij = 0. Then the equations 
become homogeneous in \l/i,2; and so their nonzero solutions exist only when 

1 - 621^12 = 0. (4) 

To be precise, for the parameters bij that do not obey the Eq. (^), the envelopes \l/i,2(0) 
are as small as r^'(O), i.e., ^1^1,2(0) ~ a/A ^ 0; this is the accuracy, under which the simple 
boundary conditions ^1,2(0) = are applicable. They are not applicable if Eq. (^) holds. It 
can be shown that the condition means the proximity of a certain interface level to the 
band edge. The energy position of this level is determined by the parameters bij and fij, 
i.e., by structure of the interface. 

Thus, the simple boundary conditions \l/i,2(0) = can be used at a sharp interface in 
the absence of interface levels close to the band edge. Otherwise, the general boundary 
conditions (|) should be used. 

It should be noted that assumption of the large bands offset at the interface is not 
important for our consideration. The boundary conditions (|I]) could be used in that case. 
This means that the simple boundary conditions \I'i 2 = can be used at a sharp interface 
even in the absence of real potential barrier there, provided that the interface levels are not 
close to the band edge. In that case the quantum confinement arises because the resonant 
tunneling of electrons is no longer possible through the interface.^ 

It is possible to rewrite the Eqs. in the more simple form: 
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til ti2 \ / ^2 
t21 t22 



(5) 



where 



_ 612621T2I - T^ll _ &12&2in2T"21 - riif22 
^11 — Z Z 5 ^12 — ~ Z Z ? 

b2i{rii - T21) &2i(tii - T21) 

. _ 1 - &12&21 , _ buhiTu - T22 

t2i — ~ 2 Z 5 ''22 — ~ Z Z • 

&2i(t-11 - 7-21) fe2l(ni - T21) 

Then the Eq.(D takes the form ^21 = 0. 

The interface influence on the electrons is determined by the parameters tij. They are 
not independent. It follows from the Eq. (|^) that the determinant of the | \tij\ \ matrix is equal 
to mi/m2. Moreover, tu ~ a, and so it is possible to assume ti2 = by the appropriate 
choice of the position of the plane z = within the unit cell at the interface. 

Thus, there are two parameters, tu ~ 1 and ^21 ~ cl~^, that determine the interface 
influence on the electron. One of them, ^21? is sensitive to the position of the interface level: 
it vanishes when this level coincides with the band edge. Another one, tu, can be considered 
as a trial parameter. 

To consider the hole states in the spherical quantum dot, we write the Luttinger Hamil- 
tonian in the spherical approximation as follows:0 

H=U + ll)^-^{pjf, (6) 
V 2 J IniQ rriQ 

where p is the momentum and J are the 4x4 matrices of the angular moment J=3/2; 7 > 
and 7i are the Luttinger parameters that relevant to the light and heavy effective masses of 
the holes: mi = mo(7i + 27)"-'^ and rrih = "^0(71 ~ 27)"^, mo is mass of the free electron. 

The moment F = 1/2, 3/2, and its projection M are the good quantum numbers due 
to the sphericalsymmetry. Solutions of the Schrodinger equation with the Hamiltonian 
are of the formEJ 



X 



mil 

(I 3/2 F \ 
M j Wigner symbols, and Xi^ is the eigenvector of the matrix. 

The radial functions Rf,f+i/2 and Rf,f-3/2 that relevant to the even solutions obey the 
equation^ll 

(71 - 2'y cos aF)PFPFRF,F+l/2 + 2^ sin aFPFP-FRF,F-3/2 + ~^\.^ ~ U{r)]RF,F+l/2 = 0, 
(71 + 27 COS aF)P^FP-FRF,F-3/2 + 2-f siu aFP^FPFRF,F+l/2 + '^[E - U {r)]RF,F-3/2 = 0, (7) 



2F-3 . ^ „ d F + 3/2 „ , d F -1/2 

— sm«p>0, PF = -r + P^ = -, '-. 

Ap dr r dr r 
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Let us, at first, suppose that Eqs. (J^) hold also at the interface where the potential 
U{r) restricts the hole movement. Then we can obtain the boundary conditions for the 
radial wave functions. Two of them arise after integration of the Eqs. (|^) over the narrow 
region |r — ro| < w/2 (a ^ w ^ A) at the interface. To obtain another two boundary- 
conditions, we have to multiply the Eqs. (|^ by r— ro before the integration. After elimination 
of RF,F+i/2{ro + w/2) oc exp{--fF+i/2r) and RF,F-3/2{ro + w/2) oc exp( -7^-3/2?") (where 
'yF+1/2 > and 7F-3/2 > are the decay exponents of the wave functions off the dot 
boundary) from the derived equations, we obtain 

{-fF+i/2W+ + V+){-fi + 2-fcosaF) - 27sin Af(7f+i/2W^- + V_) 

-lF+l/2riF,F+l/2 + rCpF+i/2 - ^2 _ 4^2 ' 

_ {-fF-3/2W- + V-){-fi - 27 COS Of) - 27 sin Xf{jf-3/2W+ + V+) 

-7-F-3/2^F,F-3/2 + ^F,F-3/2 " ^2 _ 4^2 • 

Where V+ = 2mo/ri/2t^(r - ro)RF,F+i/2dr, V- = 2moJ^l^/^U{r - ro)RF,F-3/2dr, W+ = 

-2mo j'^i^/2{r - ro)t/(r - ro)i?F,F+i/2 dr, and W- = -2mo /ri/2(^ ~ ro)U{r - ro)/?F,F-3/2 dr. 
These values vanish when w — 0, if the potential U{r) has not any singularity at 
the interface. This leads to the simple boundary conditions Rf,f+i/2{—^/1f+i/2) = 

RF,F-3/2{-fl/lF-3/2) = 0. 

The values of V± and W± don't vanish at a sharp interface where the potential U{r) 
changes essentially on the scale of the lattice constant. In particular, this is possible at a 
strain interface due to mismatch of the lattice constants of the bordering materials (e.g., at 
the Ge/Si interface). Then U{r) can be estimated as ?7 ~ D{6a/a) for \r — ro\ < a, where 
D ~ 10 eV is the constant of the deformational potential and Sa is the lattice mismatch. So 
that W ^ 6a/a ^ 1 and V ~ 6a/ a^. The boundary conditions (|]) accept the form of the 
Eqs. (^ after expansion of the radial wave functions in the integrands. This is the case even 
when 7F-3/2 and 7f+i/2 oo. 

We shall use the boundary conditions (0) in the form (^ to obtain the hole spectrum of 
the quantum dot. The radial wave functions in the free space (U = 0) arei 

RF,F+i/2ir) = AjF+i/2ikr) + BjF+i/2{kr\[^), (9) 

RF,F-3/2{r) = AijF-.3/2ikr) + BijF-3/2{kr\fp), 

where ii{z) are the spherical Bessel functions, Ai = Ataii{aF/2), Bi = i? cot(ai?/2), 
costtF = {2F — 3)/(4F), sin cti? > 0, /5 = mi/rrih] A and B are the constants that are 
determined by the boundary conditions at r = tq. 

By substitution (|]) into the boundary conditions (|^) we obtain the system of equations 
which is homogeneous in A and B. Its nonzero solutions exist only when the determinant 
vanishes, i.e., 

hijF+i/2{.krQ) - tan ^jF-3/2{kro) t2ijF+i/2{kro^J]3) + t22fF+i/2i.kro\fp) + cot '^j'F-3/2i.krQ\[p) 
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t2ijF+i/2{kro) + t22j'F+i/2{kro) - tan ^j'F-3/2{kro) tii3F+i/2{krQ\[p) + cot ^jF~3/2{kro^) 
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Where j' = dj /dr. The Eq. ([T0| ) determines the hole spectrum of the quantum dot: En = 
(71 — 2'y)h'^k'^/2mQ, where kn are the roots of Eq. ([To|) . Influence of the interface on this 
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spectrum is determined by the parameters tij. To estimate the energy Eq of the interface 
hole state, we assume k = m, where ktq ^ 1. Then from the Eq. ([To[) we obtain 



t2i (tan ^ + cot 


) 


(tnt22 - 1)(1 - v^) + til (v^cot ^ + tan ^) 


-t22 1 


(v^tan^ + cot ^) 



(11) 



So that Eq = —h'^K?/2mh. The simple case that corresponds to \I'i,2(0) = follows from the 
Eq. ( pup if we assume there ^21 00. This is possible when ^21 ^ k. The value of ^21 can 
be estimated from Eq. (|ll|), ^21 ~ ^ = ^~^y^2m/i|£'o|- Therefore the boundary conditions 
^1,2(0) = are applicable at a sharp interface, if \Eq\ ^ h^k'^ /2mh, i.e., when the energy 
of the interface level much exceeds the energy of the hole. Otherwise, the general boundary 
conditions (|) should be used. 

Figure 1 displays the left side of the Eq. (|TUp as a function of kr^. We assume tn = 1, 
m/j = mo, P = 0.1 and obtain ^21 from the Eq. (|10D provided Eq = 0.01 eV. The dashed 
curve presents the similar dependence that follows from the simple boundary conditions 
^1,2(0) = O.i We reveal an essential difference between the hole spectra. Apart from an 
essential change of the position of the roots of the Eq. ([lOD , we find that some of them 
become complex [kr^ = 15.5 ±1. Si and kvQ = 25.3 ±0.7^ on Fig. 1), and so the relevant hole 
states become quasistationary and bounded at the interface. This could be essential for the 
optical properties of the quantum dot. Moreover, such states affect the electron transport 
in the array of the quantum dots; they increase the effective cross section of the quantum 
dot. Note that the solid curve becomes close to the dashed one when Eq is about a few eV. 

The hole spectrum was found to be sensitive to the energy position of the interface level; 
namely, whether or not it is close to the band edge. Such levels really exist at the top of 
the valence band in some semiconductor /insulator contacts;lll they are responsible for the 
Fermi-level pinning. It seems that the electron interface level should be close to the valence 
band at least in wide-gap semiconductors. The interface level becomes empty then it is 
shifted too far off the top of the valence band. This results in a large surface charge and a 
strong band bending that is not favorable from the energetical point of view. Nevertheless, 
the interface level can be shifted as the result of the structure reconstruction of the interface. 
Such reconstruction does not essentially affect the interatomic spaces or angles, but it makes 
the interface level to be closer to the top of the valence band. 

In conclusion, we propose the general boundary conditions for the envelope wave func- 
tions to investigate the hole spectrum of the spherical quantum dot. We show that usually 
applicable boundary conditions ^1^1,2(0) = can be used at a smooth interface or at a sharp 
one provided that the energy separation of nearest to the band edge interface level much 
exceeds the energy of the hole under consideration. Two real parameters are sufficient to 
determine an interface influence on the hole spectrum. They could be measured in optical 
experiments or estimated theoretically [e.g., from Eq. (H)] for a certain model of the interface 
structure. The boundary conditions (|^) can be used also to describe the intervalley mixing 
of the electron in the conduction band. 
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FIGURES 




FIG. 1. Left side of Eq. ( p!o[) as a function of kvQ (bold curve). Similar dependence which 
relevant to the boundary conditions ^'1,2(0) = (dashed curve). 
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